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ABSTRACT 

The imaging channel on the Mid-Infrared Instrument (MIRI) is equipped with 
four coronagraphs that provide high contrast imaging capabilities for studying 
faint point sources and extended emission that would otherwise be overwhelmed 
by a bright point-source in its vicinity. Such bright sources might include stars 
that are orbited by exoplanets and circumstellar material, mass-loss envelopes 
around post-main-sequence stars, the near-nuclear environments in active galax¬ 
ies, and the host galaxies of distant quasars. This paper describes the coron- 
agraphic observing modes of MIRI, as well as performance estimates based on 
measurements of the MIRI flight model during cryo-vacuum testing. A brief 
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outline of coronagraphic operations is also provided. Finally, simulated MIRI 
coronagraphic observations of a few astronomical targets are presented for illus¬ 
tration. 

Subject headings: High Contrast Imaging, Astronomical Instrumentation, Extra¬ 
solar Planets, Circumstellar Disks, Quasar Host Galaxies 


Introduction 


High contrast imaging has made major strides since the first coronagraphs were used by 


Lyot to investigate the Solar Corona (Lyot 1939). The success of coronagraphy using the 


Hubble Space Telescope showed the clear advantage of the very stable environment 

achieved with space-based platforms, while the complementary successes of ground-based 
systems showed the importance of achieving very small Inner Working Angles (IWAs). 

The James Webb Space Telescope (JWST) was recognized during its early development 
as a powerful platform for studying exoplanets, or more generally, circumstellar environ¬ 
ments. In consequence, coronagraphic devices were foreseen as integral components of the 
science instruments. The potential is particularly large in the deep thermal regime where 
the Mid-Infrared Instrument (MIRT) operates, wavelengths where very faint objects can be 
detected by JWST in comparison with the limited sensitivity of even the largest ground- 
based telescopes. Therefore, in addition to a classical Lyot coronagraph, MIRI incorporates 
one of the most innovative concepts, the four-quadrant phase mask coronagraph (4QPM, 


Rouan et ah 2000) to achieve this capability. At the time of the MIRI conceptual design 


phase, only a few coronagraphic technologies had been validated both in the laboratory and 
on-sky (Baudoz et ah 2006 Boccaletti et ah 20041); the 4QPM appeared to be the best 


option to provide the smallest possible IWA of 

10 - 16 fim). 


IX/D at MIRI wavelengths (in this case. 


Direct imaging of exoplanets is one of the most challenging topics in modern astrophysics 
as it requires a large instrumental contrast, and hence very good optical quality. By the time 
of JWST’s scheduled launch in late 2018, SPHERE (Spectro-Polarimetric High-contrast 


Exoplanet Research, Beuzit et ah 2008) and GPI (Gemini Planet Imager, Macintosh et ah 


2008) will have been operating for nearly hve years. These near-IR (0.95-2.3/im), eXtreme 


Adaptive Optics (XAO) facilities will have likely discovered dozens of young giant planets 


Mhe Near Infrared Camera and Multi-Object Spectrometer (NICMOS), Space Telescope Imaging Spec¬ 
trograph (STIS), and the Advanced Camera for Surveys (ACS) all provided coronagraphic capabilities. 
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that JWST will be able to observe at complementary wavelengths (5-28/im), improving the 
spectral characterization of the planet atmospheres and measuring their energy balance. In 
addition, the Large Binocular Telescope Interferometer (LBTI) operating in the thermal 


infrared (Hinz 2009) will provide insights to the behavior of planetary systems at higher 


angular resolution than can be achieved with MIRI on JWST. 


In addition, MIRI coronagraphy will enable a range of investigations to measure struc¬ 
tures in circumstellar disks, potentially identifying locations of forming planets, or the in¬ 
fluence of unseen planets on the circumstellar material. It will also enable detailed studies 
of the circumnuclear environments in active galaxies, potentially probing the connection be¬ 
tween star-formation and outflows powered by the central supermassive black holes. Thus, 
high contrast capabilities provided by MIRI in the mid-infrared will open new frontiers in a 
variety of astrophysical contexts that have never been accessible before. 


This paper updates the descriptions of the MIRI coronagraph performance in Boccaletti 


et ah (2005) and Cavarroc et ah (2008). Section introduces the coronagraphic modes 


available in MIRI and section describes the numerical simulations and the expected raw 
contrast performance. Section presents the contrast performance measured in the lab 
during the integration and test of the flight model. Section [^discusses the target acquisition 
procedure that is specihc for the MIRI coronagraph. The potential interest of MIRI for 
exoplanet and disk science is illustrated in section together with simulated observations of 
known targets. Finally, we indicate in sectionj^some possible expansions of the high-contrast 
capabilities of MIRI. 


2. Coronagraph Implementation 

MIRI has four individual coronagraphs, one of which is based on the classic design of 
Lyot and three of which are based on four-quadrant phase masks (4QPMs). The classical 
Lyot coronagraph places an occulting spot in the focal plane to block the light from a bright 
point source from entering the instrument. Typically, this spot is of order 3 — 6X/D in 
radius so that it blocks the majority of the light from the Airy core including a few bright 
rings. Nonetheless, some extraneous light escapes 1.) as the diffraction pattern associated 
with the telescope aperture and a particular coronagraph design, and 2.) as phase and 
amplitude aberrations of the wavefront in the optical train. Some of the unwanted light can 
be intercepted by appropriate Lyot stops, at the expense of lost light, at a reimaged pupil. 
These stops can remove the diffracted light induced both by the telescope pupil and the 
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coronagraphic occulting spolj^ However, completely suppressing the telescope diffraction 
pattern resulting from the multi-mirror JWST primary would require relatively complex 
stops. For MIRI, a Lyot coronagraph is provided at 23 /xm, selected in part because of the 
lack of suitable transmissive optical materials to fabricate a 4QPM and in part to provide a 
broad spectral band to maximize the sensitivity on planetary debris disks (currently available 
4QPMs are not achromatic). 


Although classical Lyot coronagraphs can provide excellent contrast outside the area 
blocked by the occulting spot, obviously they are limited with regard to IWA to the projected 
radius of this spot (> 3X/D). Phase mask coronagraphs are designed to reduce the IWA 
to near IX/D, by replacing the occulting spot with a transparent mask that imparts phase 
differences in different parts of the focal plane so when the reimaged pupil is formed the light 
interferes more destructively than with a Lyot mask, hence rejecting the starlight outside 
the geometrical pupil. There are several implementations of phase-mask coronagraphs based 


on this principle. The theory of 4QPMs is developed in Rouan et ah 


Rouan et al. 

... 

) and 

Rouan et 


al. (2007). 


The 4QPMs in MIRI use an optical element that retards the phase by vr in two diagonally 
opposite quadrants. If a monochromatic source is placed exactly at the center of the resulting 
four-quadrant phase mask, the rejection is formally complete. The overall operation is 
illustrated in Figure The MIRI phase masks are monochromatic versions of the 4QPM, 
that is they are effective at a single wavelength. The MIRI devices generally operate with 
spectral bandpasses of about 10%. However, the net coronagraphic rejection suffers from 
several limitations including phase and amplitude aberrations, which are larger in the case of 
MIRI than the chromaticity caused by the hlter bandwidth (Riaud et al. ||2003). As another 


defect, a point-source lying in between two quadrants is also partly attenuated (Riaud et al. 


2001). Therefore, the 4QPM reduces the sensitivity (at maximum by a factor of 10) in the 


held along the four edges of the mask (in a width of about lA/D). 


All four of the coronagraphic masks are implemented in the entrance held of the MIRI 
imager on a frame supporting four sub-helds of view; the three 4QPMs provide helds of 
24" X 24", while the free-standing Lyot spot mask is held by supporting bars across a 30" x 30" 
aperture. The mounting arrangement, along with that for the low resolution spectrograph 
slit, is shown in Figure The location of the coronagraphic helds is shown in a hat held 
image (see Figure [^; they occupy about a third of the detector square held of view. No 
mechanism is required to utilize the coronagraphs - the telescope is used to position a star 


^In a coronagraph free of optical aberrations, the distribution of the intensity in the reimaged pupil is 
determined by the convolution of the telescope pupil and the Fourier Transform of the coronagraphic mask. 
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on each one as reqnired. 


The Lyot spot mask is machined in alnmininm and has a radins of 2.16", which is 3A/D 
in radins at 23/im The 4QPMs are etched in germaninm (Fignre 4) and then anti-reflection 
coated. The pnpil masks, monnted in the hlter wheel, are shown in Figure The four 
masks are each associated with an appropriate hlter (labeled ” C” in Table together with 
an appropriate Lyot stop in the hlter wheel where the pupil image is formed. The stops have a 
transmission of 62% and 72% for the 4QPMs and the Lyot spot mask, respectively. The phase 
masks were designed to operate at 10.65, 11.40 and 15.50 /im for exoplanet science while 
the Lyot coronagraph is optimized for 23 /rm, in particular for disk imaging. The two bands 
near 11 /im were selected to allow detection of the NH 3 feature at 10.65 /im in exoplanetary 
atmospheres. The choice of hlter wavelengths for MIRI coronagraphs is discussed further 


m 


Boccaletti et ah (2005). In practice, the central wavelengths of manufactured hlters 


measured at nominal temperature are well controlled ( 1 % ohset at maximum), while the 
bandwidths can be broader by 40% at maximum. 


A summary of coronagraphic conhgurations is given in Table The rejection is calcu¬ 
lated on-axis and accounts for all the aberrations considered in the simulation (see section 
1^. The coronagraphs produce a radial attenuation plotted in Figure which quantihes 
the total attenuation of a point-source located on-axis. The 50% transmission radius, de- 
hned as the IWA, is about lA/D for the 4QPMs and 3.3A/D for the Lyot spot mask. 100% 
transmission (relative) is achieved at 3 and 5A/D for respectively the 4QPM and the Lyot 
spot. 


3. Laboratory tests 
3.1. Prototype units 


4QPM coronagraphs were tested in the laboratory by 

Riaud et al. 

I2(i(i;; 

) and demon- 

strated on the telescope ( 

Boccaletti et al. [2004 

2009 

2012 

Gratadour et al. ||2005 

Riaud et 


al. 2006). The laboratory results showed very high performance (raw contrast approaching 


10 °), although use on ground-based telescopes was challenging because of the level of aber¬ 
rations left uncorrected by adaptive optics as well as the need to maintain the star precisely 
on the center of the phase mask, a near-corollary of the small IWA. 

This issue will be mitigated by operation in space with JWST. Therefore, the focus has 


shifted to tests of prototypes of the 4QPM coronagraphs for MIRI as described by [Baudoz 
et al. (2006) and Cavarroc et al. (2008). The hrst of these papers concentrates on 


the specihcations and manufacture of the 4QPM itself, optimized for the mid-infrared as 
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required for MIRI. It considers the requirements on the chromaticity of the mask material, 
how well its dimensions must be matched to the wavelength of operation, and the hdelity 
that must be achieved at the intersection of the four quadrants (where the star image must 
be placed). It shows the results of a variety of manufacturing possibilities and demonstrates 
good results using germanium as the mask material. The steps between quadrants on this 
successful mask were reactive ion etched. The second paper is concerned with the operation 
of 4QPMs specihcally in MIRI. It describes measurements of the performance of pre-flight 
phase masks in the MIRI imager test system at CEA-Sacley, including practical effects such 
as a mis-alignment of the instrument pupil with the telescope exit pupil. These experiments 
led to further optimization of the coronagraph. The performance with the flight units (next 
section) was predicted accurately in this work. 


3.2. Verification in flight instrument 


The performance of the flight coronagraphs was verihed in the flight model imager prior 
to its delivery to be integrated with the full MIRI Optics Module (OM) (|Ronayette et ah 


2010). The test arrangement of a warm source outside the cryostat with heavy Altering 


to block room-temperature thermal radiation is described in Paper III. A 30 pm diameter 
pinhole source (much smaller than the FWHM of the diAraction pattern) in the external 
apparatus was scanned in a raster pattern to measure the response of the coronagraphs, 
with results shown in Figure The faint extended residual in the image to the right is due 
to local heating of the surroundings of the pinhole infrared source, which is not rejected by 
the coronagraph and places an artiflcial floor on its contrast performance. 


The coronagraph shows good rejection of the central point spread function (PSF), by 
a factor larger than 100 on-axis (up to a few hundreds), although the level of rejection is 
probably limited more by the experimental limitations than by the coronagraph. The scans 
also showed that the IWA, defined as the angle with 50% of the maximum transmission, 
is at 1.3 A/D for the 4QPM units. Figure shows the measured contrast as a function of 
distance from the center of the PSF, and demonstrates that it agrees closely with theoretical 
expectations. 
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Simulation and Performance 


4.1. Coronagraph contrast performance 

With the performance of the coronagraph itself conhrmed, we now turn to simulations of 
its operation in JWST. Two models are under evaluation: Case A makes relatively optimistic 
assumptions about the wave front error, while Case B is more conservative. 

The simulation of MIRI coronagraphy under Case A is performed in two steps, hrst 
the calculation of coronagraphic images and PSFs themselves and then a photometric treat¬ 
ment that includes sources of noise. The images are modeled with a Fourier-based approach 

(2005)]^ The model accounts for the telescope pupil (18 


described in Boccaletti et ah 


hexagonal segments separated by small gaps) and phase aberrations (independent on each 
segment), for a total of ~130nm root mean square (rms) (mid and high frequencies) wave- 
front error (WFE), as well as a global low order aberration in the form of defocus equivalent 
to 2mm at the telescope focal plane. In the pupil plane (hlter wheel), the position of the 
Lyot stop can be mis-aligned so it is offset by 3% of the size of the pupil and includes 
also a 0.5° rotation with respect to the perfect pupil orientation. To model the telescope 
pointing variations the calculation of the coronagraphic image is repeated 200 times with 
a 1-sigma dispersion of 7mas (corresponding to the current specihcations for JWST) and 
all realizations are summed in the end. A reference star image is calculated simultaneously 
(same setup as the target) but considering different jitter realizations. We assume that the 
target/reference duty cycle can be short enough to neglect the variation of the wavefront 
errors and consider that the target to reference mismatch is driven by the pointing precision 
of 5mas (the JWST pointing requirement at 1-sigma/axis). Finally, several off-axis images 
(PSFs) are generated in the held at various positions (0.1", 0.2", 0.5", 1", 10") to account for 
the coronagraphic radial attenuation (see Figure]^. Both types of coronagraphs are included 
in the simulation, the 4QPM and the Lyot spot mask. For the broad band hlter F2300C, 
polychromatic images are generated according to the tabulated transmission of this hlter; the 
4QPM coronagraphs were simulated with their integral hlters. Raw coronagraphic contrasts 
(azimuthally averaged), corresponding to the hrst step of the simulation (free of detection 
noise), are shown in Figure]^ together with 5-a post-processing contrasts after reference star 
subtraction. This last step provides an estimate of the maximum contrast achievable directly 
with this technique. However, in addition to this simple, single reference star subtraction. 


several more advanced post-processing techniques have been developed (Lafreniere et al 


2007 Soummer et al. 2012) and should be implemented to improve this limit of detection. 


^The simulation software was used formerly to derive specifications of the MIRI coronagraphs and to 
estimate performance. 












We now consider Case B, which assnmes a WFE of 204 nm and the specihcation offset¬ 
ting accnracy (dispersion of 7 mas 1-a per axis with different realization for the star and the 
reference). Jitter is not inclnded. Model resnlts are shown in Figure]^ While the achieved 
contrast is a fairly strong fnnction of optimal positioning within the coronagraph, it is even 
more sensitive to the relative locations of the science target and the reference star (see, e.g, 
Fig. 1: Cavarroc et al. 2008). The comparison with the Case A simnlation reflects onr 


current uncertainty about the exact parameters that are applicable and how their effects 
change with wavelength (e.g., the current work indicates a greater difference between Cases 
A and B at 15.5 /im than at 11.4 /rm). Work is continuing to understand these differences 
and to use the knowledge gained to optimize coronagraphic observations. 


4.2. Full system performance 


The second part of the simulation includes the detection noise. We illustrate the proce¬ 
dure for Case A. The photon noise in the coronagraphic image depends on the star distance 
and spectral type, telescope area (25m^), telescope and optics transmission (85%), the detec¬ 
tor spectral response (i.e., quantum efficiency maximum of ~80% at ~15/rm ), the Lyot stop 
transmission (62% for 4QPM, 72% for Lyot spot and 100% for Lyot bar), and the integration 
time. The net transmission in the coronagraphic filters is about 15-20%. The intensity of the 
most distant off-axis PSF at 10" serves as a normalization factor for all simulated images, as 
if it were an observation of the unocculted star. In addition, we assume that the reference 
star is identical to the target star in terms of magnitude and spectrum. 


At the level of the detector, we included read out noise (20e- RMS), and flat field sta¬ 
bility (0.1%). The thermal background accounts for several components: the zodiacal light 
dominates at A < 15/rm, while at longer wavelengths the telescope (primary and secondary 
mirrors) as well as the deployable tower assembly shield and sunshield are mostly respon¬ 
sible for the thermal emission. Overall the model corresponds well to the specihcations of 
3.9MJy/sr at lO/rm and 200MJy/sr at 20/im (the possible higher emission from the sun- 
shield than originally expected will compromise the performance of the Lyot system, but has 
less effect on the shorter wavelength coronagraphs). 


Figure displays four cases for various filters. Each sub-panel shows the averaged 
contrast for the PSF (blue line) and the coronagraphic image (red solid line) as well as the 
3-cr contrast for the reference star subtraction process (red dashed line) and the ideal noise 
level if only photon noise dominates (red dash-dotted line). Typical contrasts after post¬ 
processing achieve 10“^ to 10“® for separations larger than 0.5" — 1". Note that the radial 
transmissions of the MIRI coronagraphs are not accounted for in these contrast profiles. 






9 


The performance presented in this section is strongly related to the end-to-end optical 
system assumptions, such as the wavefront error and the whole optical train stability of the 
telescope. The simulated performance may be conservative, particularly with the develop¬ 
ment of operations and analysis techniques optimized for JWST. In any case, the current 
model will have to be updated once the telescope is on orbit and commissioned. 


5. Operations and Target Acqnisition 


MIRI operations, including coronagraphy, are discussed in Gordon et ah (2014: Paper 
X). Here we briefly discuss the operations unique for coronagraphy, particularly target ac¬ 
quisition (TA). To achieve maximum contrast, all coronagraph observations will require TA 
to center the target accurately on the coronagraph. The 4QPMs are particularly sensitive 
to slight offsets from the optimum location (Figure 1: Cavarroc et ah 2008). The Lyot 


coronagraph is less sensitive to small offsets, but still requires precise, and repeatable TAs. 


5.1. Target acquisition 

The basic approach to coronagraphic TA starts by selecting a suitable hlter. For faint 
objects, the hlter can simply be that associated with the coronagraphic mask, but for very 
bright objects, a neutral density hlter will be needed. The target is initially placed at a 
hducial location within the coronagraphic held of view. An exposure is obtained, a centroid 
is found for the target, and the ohset necessary to move the target to the optimal location 
at the center of the coronagraph is calculated. The observatory then makes a small angle 
maneuver to place the target at the proper location in the coronagraph. Once the maneuver 
is complete the observatory enters “hne lock,” which maintains precise pointing control 
throughout the observation. If the neutral density hlter was used for TA, then the correct 
mask hlter is selected for observations (the ordering of this step is under study, in case it is 
advantageous to move to the mask hlter before maneuvering to the coronagraph position). 
A TA must be executed for observations obtained with each coronagraph separately. Thus, 
a science program that requires observations using multiple coronagraphs will need to take 
these additional TAs into account. 

Two ehects make the TA process complex: 1.) for the 4QPM coronagraphs, the phase 
mask can distort the image of a star close to its center and undermine the accuracy of the 
centroid determination; and 2.) the detector arrays have latent images (Paper VII) that 
could mimic planets or other exciting astronomical phenomena if the centroiding process left 
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them close to the target star. These effects would make adequate TA very difficult at the 
nominal JWST offsetting accuracy specifications. Fortunately, it is projected that small angle 
offsets up to 20" distance will be much more accurate than dictated by the requirements. 
For example, an offset of up to 20" is expected to be accurate to 5 mas (l-a per axis). 
Simulations of the centering accuracy on the coronagraph using the projected performance 
and a fiducial distance of 2" from the coronagraph center indicate a scatter of ~7 mas (rms) 
and average centering errors of 2 - 4 mas (Lajoie et al. 2014). The details depend on the 


particular strategy, i.e., whether one utilizes a single position for target acquisition, or uses 
more than one to acquire additional information about the pointing. None of the strategies 
quite reaches the desired centering performance for the 4QPM coronagraphs (the Lyot is 
much more relaxed in this area), so further optimization is in progress. 


5.2. Coronagraphic observations and results 

After the TA process is complete, images are acquired until the necessary on-source 
integration time is achieved. The details of the actual operation of the detectors, their 
calibration and data reduction are described in Ressler et al. (2014, Paper VIII) and Gordon 
et al. (2014, Paper X). 

To maximize contrast, most coronagraphic observations will require coronagraphic ob¬ 
servations of a “reference star,” whose image is then subtracted from the science target 
image(s). Typically this will be a star that is of similar or greater brightness compared 
with the science target; this ensures equal or greater signal to noise for the reference star 
in a similar or shorter amount of total integration time. The procedure for obtaining these 
observations is identical to that for the science target, including the TA. 

As mentioned in Section 4, several techniques and algorithms for improving contrast 
have been developed that are more sophisticated than the classical reference star subtraction 
method. One is angular differential imaging (ADI), where the target and reference star are 
observed at multiple orientations of the observatory. Because the diffraction patterns and 
detector artifacts are fixed to the telescope reference frame, they move relative to the target 
when the observatory is rolled. Subtracting two (or more) rolled images removes, or at least 
reduces, these diffraction and detector features. JWST has limited roll capability (±5°), 
because it must keep the sunshield correctly positioned relative to the Sun, but will still 
provide sufficient angular roll for some science cases at large enough angular separations 
from the star. Larger roll angles can be achieved by observing the objects at different epochs 
during the year albeit with more significant changes in the wavefront error compared to 
back-to-back observations. Another class of enhancement is centered on sophisticated image 
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analysis techniques, such as LOCI or KLIP (Lafreniere et ah 2007; Soummer et al. 2012). 


6. Illustration of Dedicated Science Cases 


This section presents a few science cases making use of the capabilities of the MIRI 
coronagraphs. We compare the simulated performance metrics presented in Section with 
the expected flux densities from giant planets orbiting solar type stars out to 10 pc. We also 
present simulated observations of the emblematic HR 8799 planetary system discovered by 


Marois et al. (2008), and observations of the debris disk associated with HD 181327. These 


examples illustrate the promise of the MIRI coronagraph in the context of exoplanet and disk 
science. In the case of exoplanets, MIRI coronagraphic imaging will: 1.) measure effective 
temperatures; 2.) determine bolometric luminosities; 3.) measure the ammonia absorption 
feature (sensitive to the presence of ammonia clouds and to the atmospheric temperature); 
and 4.) in combination with shorter wavelength measurements, support comprehensive mod¬ 
eling of the atmospheric properties (e.g., Bonnefoy et al. 2013). In the case of debris disks, 
the capability will let us: 1.) explore the virtually unknown region between about 1 and 
30 AU, containing for example the ice lines (not resolved with Spitzer, too low in surface 
brightness for HST or groundbased telescopes); 2.) constrain the minerology of planetary 
debris; and 3.) search for structures maintained by larger bodies in planetary systems. 


6.1. Exoplanet imaging 


Colored symbols in Figure pR] stand for the intensity of Jovian planets orbiting G and 
M stars at 10 pc and for different assumed planetary temperatures. The radial transmission 
of the coronagraph for each modeled planet is highlighted with dotted lines. 


For synthetic photometry of planets, we used the BDcond2000 models of [Allard et al. 


(2001) integrated over the MIRI Alters. Several test cases were considered: 300 K, 400 K, 


500 K with a radius of 1.0 Rj and 500 K, 1000 K, 1500 K with a radius of 1.4 Rj. 

At a separation of 1", the contrasts expected, even for simple reference star-subtracted 
coronagraphy, enable detection in F1140C of 1.0 Rj planets with temperature of 500 K (re¬ 
spectively 350 K) orbiting a GOV (respectively MOV) star. The situation is more favorable 
in F1550C where a 400 K temperature is reached at 1" from a GOV star. For example, the 
Pic b planet ( Lagrange et al. [ 2010) has a temperature of 1600-1700 K and a large radius 
(1.3-1.6Rj, Bonnefoy et al. 2013), hence should be brighter than any planet in Figure 10 


so almost visible in raw coronagraphic images (providing the debris disk itself is fainter than 
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the planet). Post-processing shonld yield a very high signal to noise detection. 


6.2. Imaging HR 8799 


For the simnlations of HR 8799, we assnmed an AOV star at 40pc, and 3 honrs of on- 
sonrce integration in each MIRI hlter. The planets b, c, and d, were assigned temperatnres 
of 900K, HOOK and HOOK with identical radii of 1.25Rj. The corresponding contrasts in 
each hlter, according to Allard et al. (2001), are listed in Table The angnlar separations 
are 1.72", 0.96", 0.62" for planet b, c and d. Some of the planets are close enongh to the 
center of the 4QPM coronagraphs to be attennated signihcantly. Planet b is always at a 
separation larger than 3A/D in any hlter so it is not impacted. Planet c is inward of 3A/D 
in F1065C, FH40C and F1550C, and hence is attennated by 10, 15 and 30%, respectively. 
Similarly, planet d is attennated by 35, 40, 40% in the same coronagraphic hlters. Simulated 


images are displayed in Figure 11 


Planets b and c are detected at all wavelengths, while planet d is only visible at 15.5p.m. 
Signal to noise ratios measured in the brightest pixel are given in Table These numbers, 
if integrated over a PSF size can be as large as 100 for planet b, for instance. The apparent 
elongation of planet d in the F1550C hlter (in which the star residuals are fainter than the 
planets and the background level) is in fact due to the presence of planet e nearby. These 
images were used to measure the planet photometry and to compare with the initial planet 
intensity. An example is shown in Figure [T^ for planet b, where the solid line represents the 
model degraded to a spectral resolution of i? = 20, red points stand for the initial photometry 
and blue points are the measurements obtained in each hlter. There is good agreement in the 
four hlters around A = 10pm. At F1550C the departure of the measured intensity from the 
true value is larger and likely the result of residual dihraction on one side and background 
noise on the other side. This does not mean MIRI will fail to obtain accurate photometry 
in these cases but rather that a special analysis will be required. 


6.3. The HD 181327 debris disk 


Another example of MIRI coronagraphic capabilities is given for the case of an ex¬ 
tended object. Figure [T^ shows simulated images of the ring-like debris disk around the star 
HD 181327 ( Schneider et al. ||2006 ). Here, we assumed an F5V central star, located at 50pc, 
and observed for 1 hour in each MIRI hlter. The disk model corresponds to a ~90 AU-wide 
belt containing 0.05 (up to 1 mm) of icy and porous, silicates and carbonaceous dust 
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grains in collisional eqnilibrinm. This model was fonnd by Lebreton et al. (2012) to be 


the best Bayesian solntion to reprodnce the disk Spectral Energy Distribntion (SED) from 
mid-infrared to millimeter wavelengths, assnming the radial brightness prohle measnred in 
newly-rednced HST/NICMOS scattered-light images. Scattering and absorption efficiencies 
as well as eqnilibrinm temperatnres were compnted as a fnnction of grain size distribntion 
and composition nsing state-of-the-art radiative transfer techniqnes. Monochromatic images 
were prodnced at the effective wavelength of each filter. They inclnde a scattering phase 
fnnction with qhg = 0.3 as seen in the HST images; their absolnte flnx is conservatively 
scaled to the HST one at 1.1 fim to correct for the limited capability to predict scattering 
properties. The transition between scattering- and thermal emission-dominated regimes oc- 
cnrs aronnd 10 /im, which tnrns ont to be the least favorable band for imaging dne to an 
nnfavorable compromise between contrast and resolntion. The corresponding star to disk 
total flnx ratio is 154, 48, 3.5, 0.16, respectively in F1065C, F1140C, F1550C, F2300C filters. 
The 86 AU ring-like pattern is obvions in all images as well as the north/sonth asymmetry. 
The 23/im image is sensitive to the onter part of the disk while the central cavity is blocked 


by the Lyot spot. Compared to the HST images obtained in the near-IR (Schneider et al. 


2006) and visible (Schneider et al. 2014), the MIRI images have greatly rednced stellar 


residnals in the center and will allow a cleaner analysis of the dnst grain properties based 
on mid IR colors. In addition, the ring disk is seen close to face-on and hence wonld be 
immnne to the gain bronght by angnlar differential imaging in gronnd-based observations. 
Here, for simplicity we assnmed that the radial transmission of 4QPMs is isotropic, so we did 
not model the cross-like pattern of this coronagraph where the field is partially obscured. In 
practice, for extended objects, observations at multiple rolls (at least two at 45°) will enable 
’’filling in” the gaps caused by the 4QPM boundaries and by the support bar in the Lyot 
coronagraph. 


7. Future Capability Expansions 
7.1. Expanding the wavelength coverage 

MIRI coronagraphy is centered on four filters, which can be regarded as a minimum set 
for spectral characterization of exoplanet atmospheres considering the large spectral range 
available. In particular, wavelengths shorter than lO/im carry the signatures of water vapor 
and methane, which are expected to be abundant in giant planet atmospheres. In principle, 
one can use the broad band filters for imaging (labeled ”W” in Table in combination with 
the Lyot mask. The absence of a Lyot stop to accompany these filters significantly reduces 
the coronagraphic rejection. However, the bar supporting the Lyot mask is sufficiently large 
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compared to the size of the Airy disk, especially at the shortest wavelengths, for some 
coronagraphic effect to occur using the bar as the focal plane occulter. At a minimum, these 
non-standard degraded coronagraphic modes should be considered as anti-saturation devices 
rather than coronagraphs, but they can still can be very benehcial for observing near bright 
stars. There are currently no plans to implement these modes on-orbit. However, we present 
them here for completeness and to advocate interest in them. 


To illustrate their application, similarly to section we performed simulations of the 
Lyot bar mask for three imaging Liters (F560W, F770W, FIOOOW). We assume the central 
source is offset by 6” from the center of the held so the star falls on the bar rather than 
the Lyot spot. After reference star subtraction, the contrast performance is comparable to 
that of the coronagraphic modes ( Figure 14) as the PSF falls off more rapidly with radius 
at wavelengths shorter than those passed by the coronagraphic Liters. However, we note 
for F560W and F770W that the coronagraphic eLect is really small beyond the IWA, as 
expected. 


7.2. Enhancing contrast with micro-offsets 


Another possibility is to obtain multiple observations of the reference star with slightly 
diLerent positional oLsets. JWST has a Lne steering mirror (FSM), which is used to com¬ 
pensate for small perturbations in pointing control. The FSM can also be used during 
coronagraphic imaging programs to perform micro-oLsets (sub-pixel), thus building up a 


library of reference star images (Soummer et ah 2014) for use in the more sophisticated 


algorithms. Use of the FSM has the advantage that it does not require any additional TAs, 
both making this mode very efficient and avoiding additional errors due to small-angle ma¬ 
neuvers. The image library can be used for PSF subtraction with a variety of algorithms 
(e.g; LOCI or KLIP, Lafreniere et al. (2007); Soummer et ah (2012)). A preliminary 


evaluation (at this stage omitting the eLects of latent images and pointing jitter, both of 
which need to be added to future simulations) indicates possible gains in contrast by up to 


an order of magnitude, when working close to the IWA (Soummer et al. 2014). This gain 


comes at the expense of the increased observing time on the reference star, but could be 
critical to success for demanding observations. 
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8. Conclusions 


The coronagraphic capability of MIRl is unique as it allows high-contrast imaging at 
wavelengths that are usually not considered for exoplanets and other high-contrast science. 
With JWST these science applications will become attractive owing to the large aperture 
and high sensitivity. The MIRI spatial resolution is comparable to those obtained in the 
near-IR in space with HST. To exploit this opportunity, MIRI includes phase mask corona- 
graphs, which are able to offer the same IWA as the NIRCAM coronagraphs but at longer 
wavelengths. Even with simple, single reference star subtraction, we have shown that the 
expected performance of the MIRI coronagraphs will open a new frontier in the characteriza¬ 
tion of many astronomical objects from exoplanets, to circumstellar disks, to the near-nuclear 
regions of nearby active galaxies and the host galaxies of distant quasars. 


For example, our performance analysis indicates that the threshold of detection in terms 
of planet temperature is very comparable to the one achieved with SPHERE, about 400K. As 
a consequence, the very same planets detected from the ground in the near IR with 8-m class 
telescopes will be potential targets for MIRI. MIRI, together with NIRCAM, will complete 
the spectral characterization and determine the energy balance of young giant planets in 
large-radius orbits. 


Since the coronagraphic suite of JWST was conceived (before 2005), several new devel¬ 
opments have occurred in the held of exoplanet direct imaging, both at the instrumental level 
but more importantly at the data processing level. The MIRI performance will be enhanced 
by taking advantage of these new methods for data reduction that have shown signihcant 
improvement with respect to older techniques (e.g., Lafreniere et ah 2007 Janson et ah 


2008 Soummer et al. 2012). 
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Table 1. Coronagraphic modes of MIRI 


Filter 

Coronagraph 

Stop transmission 
[%] 

Central wavelength 
[pm] 

Bandwidth 

[pm] 

IWA 

[arcsec] 

Rejection 

on-axis 

F1065C 

4QPM1 

62 

10.575 

0.75 

0.33 

260 

F1140C 

4QPM2 

62 

11.30 

0.8 

0.36 

285 

F1550C 

4QPM3 

62 

15.50 

0.9 

0.49 

310 

F2300C 

Lyot spot 

72 

22.75 

5.5 

2.16 

850 


Table 2. HR8799 planet’s contrast in difference of magnitude. 


planet 

F0560W 

F0770W 

FIOOOW 

F1065C 

F1140C 

F1550C 

b 

10.95 

10.85 

9.60 

9.73 

9.16 

9.02 

c, d 

10.29 

10.24 

9.11 

9.12 

8.82 

8.71 


Table 3. HR8799 planet’s signal to noise ratio. 


planet 

F1065C 

F1140C 

F1550C 

b 

17 

26 

11 

c 

8 

8 

4 

d 

2 

4 

6 


Table 4. Proposed non-standard coronagraphic modes. 


Filter 

Coronagraph 

Stop transmission 
[%] 

Central wavelength 
[pm] 

Bandwidth 

[{im] 

IWA 

[arcsec] 

Rejection 

on-axis 

F0560W 

Lyot bar 

100 

5.60 

1.15 

0.72 

210 

F0770W 

Lyot bar 

100 

7.65 

2.1 

0.72 

55 

FIOOOW 

Lyot bar 

100 

9.95 

1.9 

0.72 

50 
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4QPM complex amplitude 


Fig. 1.— Principle of operation of the 4QPM coronagraph. The telescope is represented 
by the large lens to the left, for simplicity assnmed to be illnminated nniformly and to be 
ronnd as shown by the circnlar pnpil sketch to the lower left. It forms a point sonrce image 
on the 4QPM as an Airy pattern. The phase plate modihes the wavefronts to prodnce the 
complex amplitnde pattern indicated for a sonrce centered exactly on the intersection of the 
qnadrants. The lens behind the telescope focal plane forms a pnpil where a traditional Lyot 
stop is placed to remove part of the diffracted light from the telescope pnpil throngh the 
coronagraphic mask. When the following lens reimages the focal plane onto the detector, 
the complex amplitnde pattern canses the light from the bright central sonrce to interfere 


destrnctively. Fignre adapted from Ronan et ah (2007). The layont of a Lyot coronagraph 


is identical except for the snbstitntion of an occnlting spot for the 4QPM. 
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4QPM (10.65 ^im) 4QPM(11.4|im) 4QPM (15.5 iim) 

Lyot (23 |im) 


LRS slit 


Fig. 2.— Mounting bracket for the four coronagraph image-plane masks and the LRS slit. 
The action of the spring-loaded mounting hardware is shown as an exploded view. 



30 X 30 arcsec 
Lyot Mask 


24 X 24 arcsec 
4QPMS 


- 21 - 



4.7 X 0.51 arcsec '* 74 arcsec 

LRS slit 


Fig. 3.— Flat field image: the various coronagraphic fields are located to the left (phase 
masks at the bottom, and the Lyot mask at the top), while the imager field is at the right. 













Fig. 4.— Microscopic photography of the one of the 4QPMs manufactured at CEA (Com¬ 
missariat a I’Energie Atomique, France). The central defect is about 7/im in width. As a 
comparison, the size of the JWST PSF in the focal plane at 11.40/im is AF/D=230/im. 







23 



Fig. 5.— Coronagraph pupil masks, manufactured at LESIA (the Laboratory for Space 
Studies and Asrtrophysics Instrumentation, at the Paris Observatory). The 4QPM version 
is to the left and the Lyot one to the right. 
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Separation in X/D 


Fig. 6.— Radial transmission of the MIRI coronagraphs (4QPM in red, Lyot spot in bine). 
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-0.02 0.00 0.02 0.04 0.06' ‘ 0.08- ..>O.TO 



Fig. 7.— Sample images at 11.40 /im: left without the coronagraph, right with it. The held 
of view is 8.36” (76 detector pixels) and the color scale indicates the contrast with respect 
to the PSF maximum. The level of the coronagraphic image is close to the background 
noise and the central residual is attributed to the size of the source, which is not perfectly 
point-like. 
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Seporotlon in X/D 


Fig. 8.— Contrast achieved in the test of the flight 4QPM coronagraphs (red for 10.65 /rm, 
magenta for 11.40 /rm and blue for 15.50 /rm) in the flight imager. The solid lines are for 
the unattenuated PSF and the dashed ones show the effect of the 4QPM. 
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Fig. 9.— Normalized noise-free contrasts obtained under Case A (smaller WFE, in red) and 
case B (larger WFE, in blue: only the reference star subtracted final result) in F1140C (left) 
and F1550C (right) filters on the PSF (solid) and the raw coronagraphic image (dashed). 
Estimated 5cr contrasts using reference star subtraction are also shown (dash-dotted). 
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Fig. 10.— Estimated performance of MIRI in a few selected cases, assnming Case A (smaller 
WFE) and inclnding the effects of noise. See text for details. 
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Fig. 11.— Simulated observations of the HR 8799 system in the 4QPM coronagraphic hlters, 
assuming Case A (smaller WFE). The planets are labeled b, c and d. The image is in units 
of contrast relative to the un-occulted star, as displayed by the color scale. 
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Fig. 12.— Spectral Energy distribution (expressed in contrast) of the planet HR8799 b, 


assuming Case A (smaller WFE). The solid line represents the spectral model (Allard et ah 


2001) degraded to a resolution of R = 20. The red and blue points represent the true and 


measured intensities respectively in several MIRI hlters (F560W, F770W, FIOOOW, F1065C, 
F1140C, F1550C). To get high quality measurements in the hrst three hlters may require 
implementation of the coronagraph bar observations discussed in Section 7. Horizontal bars 
indicate the width of hlters. 
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Fig. 13.— Simulated observations of the HD 181327 debris disk in the coronagraphic hlters, 
assuming Case A (smaller WFE). The image is in units of contrast relative to the un-occulted 
star, as displayed by the color scale. The cross effect from the 4QPM is not included in this 
simulation. 
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Fig. 14.— Same as Figure 


for the Lyot bar mask and assuming Case A (smaller WFE) 
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